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Abstract 

Using the form factors calculated in the three-point QCD sum rules, we calculate the new 
physics contributions to the physical observables of B c — > D*fi + u~ decay in a family non- 
universal Z' model. Under the consideration of three cases of the new physics parameters, we 
find that: (a) the Z' boson can provide large contributions to the differential decay rates; (b) 
the forward-backward asymmetry (FBA) can be increased by about 47%, 38%, and 110% at 
most in SI, S2, and extreme limit values (ELV), respectively. In addition, the zero crossing can 
be shifted in all the cases; (c) when s > 0.08, the value of Pl can be changed from —1 in the 
Standard Model (SM) to —0.5 in SI, —0.6 in S2, and in extreme limit values, respectively; (d) 
the new physics corrections to Pt will decrease the SM prediction about 25% for the cases of 
SI and S2, 100% for the case of ELV. 
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I. INTRODUCTION 



The Standard Model (SM) of interactions among elementary particles is one of the 
best verified physics theories up to now but there are many open fundamental questions 
remain unanswered within the scope of the SM. High energy physics experiments are 
designed to address the open questions through the search of new physics (NP) using two 
complementary approaches. One is to discover the new particles at the high energy Large 
Hadron Collider (LHC). The other is to search for the effects of NP through measurements 
of flavor physics reactions at lower energy scales and evidence of a deviation from the SM 
prediction. 

After the observation of the rare radiative decay b — > S7 [1] , the flavor-changing neutral 
current (FCNC) transitions became more attractive and since then many works about 
rare radiative, leptonic and semileptonic decays have been intensively done in the B u ^,s 
system [2] . Among these decays, semileptonic decay channels are significant because their 
branching ratios are relatively larger. These works will be more perfect if similar studies 
for B c , observed in 1998 by CDF Collaboration [3], are also included. 

The charmed B c meson is a ground state of two heavy quarks b and c. Because of the 
two heavy quarks, the decays of the B c meson are rather different from B u j B^j B s mesons. 
Physicists therefore believe that the B c physics must be very rich compared to the other 
B mesons if the statistics reaches high level [3-5]. At LHC, around 5 x 10 10 B c events 
per year are expected [6, 7]. The expected number events are motivating to work on the 
B c phenomenology and this possibility will provide facilities to study the observables of 
rare B c decays such as branching ratios, forward-backward asymmetry and polarization 
asymmetries. 

The rare B c — > D* s l + l~ decays are proceeded by FCNC transition of b — > sl + l~, which 
are forbidden at the tree level in the SM, and play an important role in the precision test 
of the SM. Meanwhile, they offer a valuable possibility of an indirect search of NP for their 
sensitivity to the gauge structure and new contributions. Up to now, the possible new 
physics contributions to B c — > D*l + l~ decays have been studied extensively, for example, 
by using model independent effective Hamiltonian [8], in Supersymmetric models [9], with 
fourth generation effects [10], and in single universal extra dimension [11]. 

When concentrating on the exclusive B c — > D*l + l~ decays, one needs to know the form 
factors. As for B c — > D* transition, the form factors have been calculated using different 
approaches, such as light front constituent quark models [12], a relativistic constituent 
quark model [13], a relativistic quark model [14], the Ward identities [15], in light cone 
QCD [16, 17], and QCD sum rules [18-20]. In this work, we will adopt the form factors 
calculated in the three-point QCD sum rules [20] to study the Z' effects on the observables 
for B c — > D* s fi + fi~ decay. 

The general framework for non- universal Z' model has been developed in Ref. [21]. In 
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this model, Z' gauge boson could be naturally derived by adding additional U(l)' gauge 
symmetry. Non-universal Z' couplings can induce FCNC b — > s and d transitions at 
tree level. Its effects on b — > s transition have received great attention and been widely 
studied in the literature. The previous works in a family non-universal Z' model boson 
redound to resolve many puzzles, such as "ttK puzzle" [22, 23], anomalous B s — B s mixing 
phase [24, 25] and mismatch in Afb(B — > K*fi + fi~) spectrum at low q 2 region [26, 27]. 
Motivated by this, we will study the effects of the Z' boson on the rare decay B c — >■ 
. 

This paper is organized as follows. In Section 2, we present the effective Hamiltonian 
responsible for the b — > sl + l~ transition in both the SM and the family non-universal Z' 
model. In this section we also present the matrix element, and the expressions of various 
physical observables in Z' model. In Section 3, we show the numerical results of the 
observables for B c — > D*fi + fi~ decay in the SM and Z 1 model. The final section is the 
summary. 



II. EFFECTIVE HAMILTONIAN, MATRIX ELEMENTS AND OBSERVABLES 

FOR b -»• sl+r DECAY 

At quark level, the rare semileptonic decay b — > sl + l~ can be described in terms of the 
effective Hamiltonian which is given by [28, 29] 



AG, 



* i=i 



(1) 



here the explicit expressions of Oj can be found in Ref. [28], in which 

2 2 

o 9 = - 2 {d ltl p L b){m), o 10 = ^(d^p L b)(h^ 5 i). (2) 

9 S 91 

The Wilson coefficients Cj can be expanded perturbatively [30-33]. The effective coeffi- 
cients C7 9 , can be written as [28] 

reS _^ r l r 4 20 80 
C 7 -- C 7-3C 3 --C 4 --C 5 --C 6 , 

Air 

Cf = —C 9 + Y(s), 
a s 

C10 — — C10 , (3) 
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where the perturbative part Y(q 2 ) stands for the matrix element of four-quark operators 
and is given by 

Y(q 2 ) = h(m c ,s)(lc 1 + C 2 + 6C3 + 6OC5) 

1 4 64 

+-h(l, s) ( - 7C 3 - -C 4 - 76C 5 - -C 6 ) 

1 4 64 

+-/i(0, s) ( - C 3 - -C 4 - I6C5 - — C 6 ) 

4^ 64^ 64^ /A . 
+ 3 C3 + Y C5 + 27 a - (4) 

with s = q 2 /fn 2 Bc , rh c = m c /m Bc . We have neglected the resonance contribution. For the 
detailed discussion of such resonance effects, we refer to Refs. [9-11]. 

Exclusive decay B c — > D*fi + fi" is described in terms of matrix elements of the quark 
operators in the effective Hamiltonian over meson states, which can be parameterized in 
terms of form factors. The matrix elements of B c — > D* transition are given by [34] 

(D* s (p)\(V-A),\B c (p Bc )) = -ie*(m Bc + m Dt )A (s)+i( PBc+P ) ll (e* PBc ) A+ } s) 

m B c + m D* 

+iq^*PB c ) 2 ^A_{s) + W .e>^ 2Av } S) . (5) 



and 



(£>:(p)|sa^g"(l + 75 )6|B c (pB c )) = ie^ v f B f 27\(s) 

+ T 2(s) {e;(m| c - m 2 Dt ) - (e*p Bc ) (p Bc +p)^} 



Tz{s){e*p Bc ) { g M - 2 S 2 (p Bc + p)„ J> (6) 

''Be "I'D* 



here s = g 2 , = (p Bc — p)^, and e M is polarization vector of the vector meson D*. 

The dilepton invariant mass spectrum for B c — > D*l + l~ decays can be expressed by [34, 

35] 

dr _ G\ a 2 m 5 Bc 



ds 2 10 7T 5 



\V*V tb \ u(s)D (7) 
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where the explicit expression of D is 



D = 



\A\ 



m 



sA(l + 2-f) + \E\ 2 s 



uls) 



4m 2 D * 



4m 2 D * 



\B\\\- 



uls 



+ 8m 2 D ,(s + 2m 2 )) + \F\ 2 (\ - 



uls) 



+ 8m 2 D *(s - Am 2 )) 



\C\ 2 {\ 



uls) 



-) + \G\ 



uls) ~ ~ 9 



+ 4mf(2 + 2m 2 D * - s) 



Re(BC*)(\ 



uls 



-)(l-m 2 



uls 



+ Re(FG*)((A - - m 2 - $) + AmfX) 



2-^-\ [Re(FH*) - Re(GH*)(l - m 2 D *)] + ^-s\\H\ 2 , 



mt 



m 



m 



(8) 



with m ; = mi/m Bc , and m D * = m D */m Bc . The kinematic variables s and w are the same 
as Ref. [34] . The auxiliary functions A , B ,C , E , F and G which are combinations of the 
effective Wilson coefficients in Eq. (3) and the form factors of B c — > D* transition can be 
found in Refs. [34, 35]. For the convenience of the reader, we present these functions in 
the Appendix A. 

The normalized forward-backward asymmetry (FBA) is defined as 



A 



FB{S 



ds 



I-i dcos9£^Sign(cos9) 



f +1 dcose-g^ 

J — 1 asdcosa 

According to this definition, the explicit expression of FBA is: 

dApB 



(9) 



ds 



-D = u(s)s[Re(BE*) + Re(AF*)} . 



The lepton polarization can be defined as: 



dS = 2 + ( L L + N N + T t) ' ] 



(10) 



(11) 



where the subscript "0" stands for the unpolarized decay case. Pl and Pt are the longi- 
tudinal and transverse polarization asymmetries in the decay plane respectively, and P/v 
is the normal polarization asymmetry in the direction perpendicular to the decay plane. 
The lepton polarization asymmetry Pj can be derived by 



dT(n = e.i)/ds — dT(h = —ii)/ds 
dT(n = ii) /ds + dF(n — — ej) / ds 



(12) 
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the results are 



s 3 3mi 



P L D = ^/i_4_L^ _ r R e (AE*) + - — 2 s Re(BF*) 



A(l - m 2 D » - s) A 2 

— ^ -Re(BG* + CF*) + -4— Re(CG*) }, (13) 

3m j}* 3m D * 



P N D = Zl^l) [lm {FG*){l + 3m% - s) 
Am D * 1 m D * 1 3 

+Im{FH*){\ - m 2 D , - s) - Im(GH*)\] 

+ 2m D *m{[Im(BE*) + 7m(AF*)] J, (14) 

n\/\m, { (1 — mh-\* — s) r _ n 

P T £> = V ^ 4sPe(AB*) + " [~Re(BF*) + (1 - m D *)Re{BG*) + si?e(5iT)] 

4Vs [ 



A 



+ — [Pe(CF*) - (1 - m^,)^(CG*) - si2e(Cff*)] ). (15) 



In the family non-universal Z' model, the flavor neutral currents arise even at tree level 
owing to non-diagonal chiral coupling matrix. Postulating that the couplings of right- 
handed quark flavors with Z' boson are diagonal, the Z' part of the effective Hamiltonian 
for b — > sl + l~ transition is described by [24] 



1G 

u% } (b -+ sin-) = -^v tb v t ; -^^(sb) v ^ v . A -^^(sb) v ^ A (ii) v+A 



V2 



B ^ B ll,-^ (Ti\ B s b B U,-„i 



' v ' v tb v ts 



+ h.c. . 
(16) 

To extract the Z' corrections to the Wilson coefficients, one can reformulate Eq. (16) 

as 



AG 

U% f {b sin') = -^V tb V* [AC 9 9 + AC' w O 10 ] + h.c. , (17) 



with 



AC'(M W ) = - 9 4^S™, 



AC' W (M W ) = 9 4^T D u R - (18) 



e 2 V t * s V tb 

2 r>L 

S ^sb 7"), 

where S{( R = (Bf( + B R ), D\ R = (Bft - B R ) with B% and Bf{ R referring to the effective 
chiral Z' couplings to quarks and leptons, respectively. The off-diagonal element B^ b 
contains a new weak phase and can be written as \B^ b \e 1 ^ . 
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TABLE I: Default values of inputs parameters used in our numerical calculations. 



m b = 4.8 GeV, m Bc = 6.28 GeV, m D * = 2.112 GeV, m M = 0.106 GeV 
\V tb V t *\ = 0.041, a = 1/137, t Bc = 0.46 x 1Q- I2 s. 



When we include the Z' contributions with the assumption of no significant RG running 
effects between M z > and M w scales, the Wilson coefficients can be written as 

Cg$(M w ) Cg$(M w ) + AC' 9A0 (M W ) . (19) 

After inclusion of the new contributions from Z' boson, the RG evolution of the Wilson 
coefficients down to low scale is exactly the same as in the SM. 



III. NUMERICAL RESULTS 



In this section, we focus on the numerical calculations of the branching ratios, forward- 
backward asymmetry and polarization asymmetries for B c — > D*fi + fi~ decay. The input 
parameters which are related to our analysis are summarized in Table I. 

For the form factors A v (s), A (s), A + (s), A-(s), Ti(s), T 2 (s) and T 3 (s), we choose 
them derived by the three-point QCD sum rules [20] , in which the parametrization of the 
form factors with respect to q 2 are as follows: 

*« = itSw (20) 

where the values of the parameters F (0), a and (3 are listed in Table II. 



TABLE II: B c -»• D* form factors in the QCD Sum Rules [20]. 



F{q 2 ) 


F(0) 


a 


a 


Ay (q 2 ) 


0.54 


-1.28 


-0.230 


A (q 2 ) 


0.30 


-0.13 


-0.180 


A+(q 2 ) 


0.36 


-0.67 


-0.066 


A-(q 2 ) 


-0.57 


-1.11 


-0.140 


Ti(q 2 ) 


0.31 


-1.28 


-0.230 


T 2 (q 2 ) 


0.33 


-0.10 


-0.097 


Ts(q 2 ) 


0.29 


-0.91 


0.007 



In the family non-universal Z' model, the Z' contributions rely on four parameters |-BjjJ, 
g , and D^p. These parameters have been constrained from the well measured decays 



TABLE III: The inputs parameters for the Z' couplings [25, 26]. 





I^Kxio- 3 ) 




^(xKT 2 ) 




SI 


1.09 ±0.22 


-72 ±7 


-2.8 ±3.9 


-6.7 ±2.6 


S2 


2.20 ±0.15 


-82 ±4 


-1.2 ± 1.4 


-2.5 ±0.9 



by many groups [24-27]. \B^ b \ and 0f have been strictly constrained by B s — B s mixing, 
B — ?> txK^ and pK decays. After taking into account constraints from B^ — > X s pp, Kpp 
and K*pp, as well as B s — > pp decays, the bounds on and are also obtained. 
For the sake of convenience, we recollect their numerical results in Table III, with SI and 
S2 corresponding to two fitting results of UTfit Collaboration for B s — B s mixing [36]. 

Recently, CDF, DO, and LHCb collaborations [37-39] have updated the CP violation 
parameter (f> s in B s system. These precise measurements will suppress the magnitude of 
b — s — Z' coupling by about 10%, and have no effect on the new weak phase 0f . However, 
the weak phase can be constrained by the data of B — > irK^ and pK decays and the 
results are consistent with the previous Refs. [25, 26]. Indeed, the quantity that is directly 
related to the decay studied here is the product of the couplings of b — s — Z' and p—p — Z', 
and the updated experimental data of B s mixing have less effect on it. According to the 
above analysis, we will adopt the inputs parameters for the Z' couplings as in Table III 
in our theoretical calculation. Meanwhile, we also choose the extreme values of SI which 
are named extreme limit values (ELV) to show the maximal effects of Z' contributions, 
and the ELV are 



\B L sh \ = 1.31 x 10- 3 , <f> L a = -79° , S™ = -6.7 x lO" 2 , D™ = -9.3 x 10" 2 



(21) 



Using the input parameters given above, we obtain the results of the branching ratios 
both in the SM and the family non-universal Z' model without resonance contributions. 

< 2.32^; 27 x 10- 7 (SM), 



Br(B c -> D* s p + p-) 



3.363£ x 10^ 7 (SI), 

2.80T^ 2 x 10- 7 (S2), 1 ] 

5.2ll°J 7 x 10- 7 (ELV). 



The theoretical errors are induced by the uncertainties of form factors. From the numerical 
results, one can see that branching ratio for decay B c — > D*p + p~ is sensitive to the Z 1 
contributions. With respect to the central value of the SM prediction, the new physics 
contributions in the family non-universal Z' model can provide an enhancement about 
45%, 21%, and 125% for the case of SI, S2, and ELV, respectively. 

Fig. 1 shows the s dependence of the differential decay rates for decay B c — > D*p + p~ 
both in the SM and the family non-universal Z 1 model using the central values of the 
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FIG. 1: The s dependence of the differential decay rates dBr(B c —> D*fi + /i~)/ds both in the 
SM and the family non-universal Z' model. The solid, dashed, dash-dotted, short-dashed lines 
show the SM prediction, the theoretical results of SI, S2, and ELV, respectively. 




FIG. 2: The FBA of decay B c — > Z)*/i + /x as a function s both in the SM and the family 
non-universal Z' model. 

input parameters. The solid line refers to the SM prediction, while the dashed, dash- 
dotted, short-dashed curves correspond to the theoretical results of SI, S2, and ELV, 
respectively. The Z' enhancements to the differential decay rate are significant in almost 
the whole region of s and strongly depend on the variation of NP parameters. 
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0.0 0.1 0.2 0.3 0.4 

FIG. 3: The longitudinal lepton polarization asymmetry of decay B c — > D*fi + fi~ as a function 
s both in the SM and the family non-universal Z' model. 

The s dependence of forward-backward asymmetry for B c — > D*fi + fi~ decay is pre- 
sented in Fig. 2. Compared to the SM results, when including the NP effects from Z' 
boson, the FBA can be increased by about 47%, 38%, and 110% at most in SI, S2, and 
ELV, respectively. It is easy to see that the zero crossing in A FB {B C — > D* s n + ji^) also ex- 
ists and Z' corrections can shift s = 0.075 in the SM to s = 0.104 in SI, and s = 0.093 
in S2, respectively. As for the case of ELV, the Z' effects on A FB (B C — > D*fi + yr) are 
more significant and can lead zero crossing to vanish. 

In Fig. 3, we plot the longitudinal lepton polarization asymmetry of decay B c — > 
D*n + fi~ as a function s both in the SM and the family non- universal Z' model. After 
inclusion of the Z' contributions, there are also apparent deviations in the values of the 
Pl{B c — > D*[i + fi~) for all the cases in Z' model from that of the SM predictions. When 
s > 0.08, the value of the longitudinal polarization asymmetry can be changed from — 1 
in the SM to —0.5 in SI, and —0.6 in S2, respectively. In the extreme case, the Z' effects 
could flip the sign of the SM predictions when s > 0.006 and the theoretical values might 
be close to zero in large momentum region. 

The transverse lepton polarization asymmetry of decay B c — > D* function s 

both in the SM and the family non-universal Z' model is given in Fig. 4. The new physics 
corrections from Z' boson are small, and will decrease the SM prediction about 25% for 
the cases of SI and S2 in low s region. However, for the case of ELV, the decrease could 
be rather large and reach 100% of the SM predictions. In addition, the sign of Pt will be 
changed in low momentum region and its values approach to zero when s > 0.037. 
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FIG. 4: The transverse lepton polarization asymmetry of decay B c — > D*fi + fi as a function s 
both in the SM and the family non-universal Z' model. 

IV. SUMMARY 

In this paper, we calculated the Z' contributions to the branching ratio, forward- 
backward asymmetry and polarization asymmetries for B c — > D*fi + fi~ decay in the family 
non-universal Z' model by employing the effective Hamiltonian with the form factors 
calculated in the three-point QCD sum rules. 

In Section 2, we presented the theoretical framework of b — > sl + l~ transition including 
the effective Hamiltonian, matrix element and the physical observables. In Section 3, we 
showed the numerical results of the observables and made phenomenological analysis for 
B c — > D*{i + fi~ decay in the SM and the family non- universal Z' model. 

As expected, the Z' contributions to the observables for B c — > D*fi + fi~ decay could 
be significant in size. From the numerical results, we found that: 

• With respect to the SM prediction, the Z' contributions to the differential decay 
rates are significant in almost the whole region of s and strongly depend on the 
variation of NP parameters. 

• The new physics enhancements to FBA could be large, and reach 47%, 38%, and 
110% at most in SI, S2, and ELV, respectively. The zero crossing could be shifted 
from s = 0.075 in the SM to s = 0.104 in SI, and s = 0.093 in S2, respectively. 
As for the case of ELV, the Z' effects could lead zero crossing to vanish. 

• The values of Pl deviated apparently from that of the SM predictions for all the 
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cases in Z' model. In high s region, the values of Pl could be changed from — 1 in 
the SM to —0.5 in SI, —0.6 in S2, and in ELV, respectively. 

• The new physics corrections to P T would decrease the SM prediction about 25% for 
the cases of SI and S2 in low s region. However, for the case of ELV, the decrease 
could be rather large and reach 100% of the SM predictions. 
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Appendix A: Auxiliary functions 

The auxiliary functions are given as follows [34, 35]: 

A{s) = —^—C e 9 ff (s)A v (s) + ^C^T^s), (23) 
1 + rriD* s 

B{s) = (1 + m D ,)C e 9 ff (s)A (s) + ^(1 - m 2 Dts )C e 7 ff T 2 (s), (24) 

c(s) = ^ct'WMS + ^rcr (r 8 (*) + , (25) 

m = —^—Ct f f A v (s), (26) 

F(s) = (l + m D *)Ct f f A (s), (27) 

G(s) = Y^—CtffA + (s), (28) 

H(S) = 2 -^C^A_(S), (29) 
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